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1.  INTRODUCTION 

A  practical  question  is  commonly 
asked:  How  long  is  an  ocean  (or 
atmospheric)  model  valid  since  being 
integrated  from  its  initial  state?  Or  what  is 
the  model  valid  prediction  period  (VPP)?  To 
answer  this  question,  uncertainty  in  ocean 
(or  atmospheric)  models  should  be 
investigated.  It  is  widely  recognized  that  the 
uncertainty  can  be  traced  back  to  three 
factors  [Lorenz,  1963;  1984]:  (a) 

measurement  errors,  (b)  model  errors  such 
as  discretization  and  uncertain  model 
parameters,  and  (c)  chaotic  dynamics.  The 
measurement  errors  cause  uncertainty  in 
initial  and/or  boundary  conditions.  The 
discretization  causes  small-scale  “subgrid” 
processes  to  be  either  discarded  or 
parameterized.  The  chaotic  dynamics  may 
trigger  a  subsequent  amplification  of  small 
errors  through  a  complex  response. 

The  three  factors  are  related  to  each 
other  causing  model  uncertainty.  Chu  [1999 
pointed  out  that  the  boundary  errors  act  as  a 
forcing  term  (stochastic  forcing)  in  the 
Lorenz  system.  This  suggests  that  the 
oceanic  (or  atmospheric)  variables  are 
sensitive  to  measurement  errors  (uncertain 
initial  and/or  boundary  conditions),  model 
errors  (discretization),  and  chaotic  (or 
stochastic)  dynamics. 

Currently,  some  time  scale  (e.g.,  e- 
folding  scale)  is  computed  from  the 
instantaneous  error  (IE)  growth  to  represent 
the  model  VPP.  The  faster  the  IE  grows,  the 
shorter  the  e-folding  scale  is,  and  in  turn  the 
shorter  the  VPP  is.  Using  IE,  evaluation  of 
deterministic  models  becomes  stability 
analysis  in  terms  of  either  the  leading 
(largest)  Lyapunov  exponent  [e.g.,  Lorenz, 
1969]  or  the  amplification  factors  calculated 
from  the  leading  singular  vectors  [e.g., 
Farrell  and  loannou,  1996  a,  b].  For 
stochastic  model,  the  statistical  analysis 
becomes  useful  [Ehrendorfer,  1994  a,  b; 


Nicolis,  1992].  The  probabilistic  properties  of 
the  prediction  error  are  described  using  the 
probability  density  function  (PDF)  satisfying 
the  Liouville  equation  or  the  Fokker-Plank 
equation.  Solving  this  equation,  the  mean 
and  variance  of  errors  can  be  obtained. 
Nicolis  [1992]  investigated  the  properties  of 
error-growth  using  a  simple  low-order  model 
(projection  of  Lorenz  system  into  most 
unstable  manifold)  with  stochastic  forcing.  A 
large  number  of  numerical  experiments  were 
performed  to  assess  the  relative  importance 
of  average  and  random  elements  in  error 
growth. 

In  fact,  the  IE  growth  rate  is  not  the  only 
factor  to  determine  VPP.  Other  factors,  such 
as  the  initial  error  and  tolerance  level  of 
prediction,  should  also  be  considered.  The 
tolerance  level  of  prediction  is  defined  as  the 
maximum  error  the  model  can  afford  (still 
keeping  the  meaningful  prediction).  For  the 
same  IE  growth  rate,  the  higher  the 
tolerance  level  (initial  error)  the  longer 
(shorter)  the  model  VPP  is.  Thus,  the  model 
VPP  should  be  defined  as  the  time  period 
when  the  model  error  first  exceeds  a  pre¬ 
determined  criterion  (i.e.,  the  tolerance  level 
e).  The  longer  the  VPP,  the  higher  the  model 
predictability  will  be.  In  this  study,  we 
develop  a  theoretical  framework  of  model 
predictability  evaluation  using  VPP,  and 
illustrate  the  usefulness  and  special  features 
of  VPP.  The  outline  of  this  paper  is  depicted 
as  follows.  Description  of  prediction  error  of 
deterministic  and  stochastic  models  is  given 
in  Section  2.  Estimate  of  VPP  is  given  in 
Section  3.  In  Section  4,  determination  of 
VPP  for  a  one-dimensional  stochastic 
dynamic  system  is  discussed.  In  Section  5, 
the  conclusions  are  presented. 


2.  Prediction  Error 
2.1.  Dynamic  Law 

Let  x(f)  =  [x<v(t),  x<2>(t),  . 


x(n>(t)]  be  the 


set  of  variables  characterizing 


dynamics  of  the  ocean  (or  atmosphere)  in  a 
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certain  level  of  description.  Let  the  dynamic 
law  be  given  by 

£  =  f(x,0  (1) 

dt 

where  f  is  a  functional.  Deterministic 
(oceanic  or  atmospheric)  prediction  is  to  find 
the  solution  of  (1)  with  an  initial  condition 
x(0  =  xo  (2) 

where  x0  is  an  initial  value  of  x. 

With  a  linear  stochastic  forcing, 
g(f)x,  Eq.(1)  becomes 
dx 

—  =  f(x,t)  +  q(t)x  (3) 

dt 

Here,  q(t)  is  assumed  to  be  a  random 
variable  with  zero  mean 

(?(*)>  =  0  (4) 

and  pulse-type  variance 

(q(t)  q(t))  =  q18(t-t)  (5) 

where  the  bracket  <  >  is  defined  as 
ensemble  mean  over  realizations  generated 
by  the  stochastic  forcing,  8  is  the  Delta 
function,  and  cf  is  the  intensity  of  the 
stochastic  forcing. 

2.2.  Prediction  Model  and  Error 
Variance 

Let  y (t)  =  [y (1)(t),  f>(t),  ...,  y<n>(t)]  be 
the  estimate  of  x(f)  using  a  prediction 
model 

-j-  =  h(y>0  ■  (6a) 

dt 

with  an  initial  condition 

y(O  =  y0  (6b) 

where  y0  is  the  initial  value  of  y. 

Difference  between  reality  (x)  and 
prediction  (y)  at  any  time  t(>t0) 
z  =  x-y 

is  defined  by  the  prediction  error  vector  and 
this  difference  at  t0 
Zo  =  x0-y0 

is  defined  by  the  initial  error  vector.  If  the 
components  [x(1>(t),  x(2>(t),  x<n)(t)]  are  not 

equally  important  in  terms  of  prediction,  the 
uncertainty  of  model  prediction  can  be 
measured  by  the  error  variance 

Var(z)=(z*Az>.  (7) 

Here,  the  superscript  't’  indicates  the 
transpose  of  the  prediction  error  vector  z.  A 
is  an  nxn  weight  matrix. 

3.  Valid  Prediction  Period 


3.1.  Indirect  Estimation  from  IE 
Growth  Rate 

Use  of  IE  is  to  investigate  the  model 
error  growth  due  to  an  initial  error, 

z(f0)  =  z0  (8) 

where  t0  is  the  initial  time.  To  evaluate 
model  predictability  becomes  to  analyze 
stability  of  the  system  to  the  given  initial 
analysis  error,  z0.  Thus,  many  instability 
theories  can  be  easily  incorporated,  such  as 
the  leading  (largest)  Lyapunov  exponent 
[e.g.,  Lorenz  1969]  or  the  amplification 
factors  calculated  from  the  leading  singular 
vectors  [e.g.,  Farrell  and  loannou,  1996  a, 
b].  Taking  a  linear  time-dependent 

dynamical  system 


dx  .  .  , 

—  =  A(r)x, 
dt 

as  an  example.  The 

exponent  is  defined  as 

,  ln((||O(M0)||) 

A  =  lim  sup - 


first 


(9) 

Lyapunov 

(10) 


where  the  propagator  $(f)  is  given  by 
[Coddington  and  Levinson,  1 955] 


t  t  r 

<S>(t, t0 )  =  I  +  J  A (s)ds  +  J  A(r)dr  j  A (s)ds  +  ... . 

t0  t0  t0 

where  I  is  the  unit  matrix.  The  only 

information  we  can  get  here  is:  the  larger  the 
value  of  X,  the  shorter  the  model  VPP. 
Usually,  the  e-folding  scale  relating  to  the  IE 
growth  rate  is  used  to  represent  VPP. 

3.2.  Direct  Calculation 
To  quantify  VPP,  we  first  define  two 
model  error  limits:  minimum  (noise  level 
(^noise)  and  maximum  (tolerance  level  e). 
The  model  prediction  is  considered 
‘accurate’  if  the  model  error  is  less  than  the 
noise  level, 

Var(z)<Use2.  (11) 

The  model  prediction  is  meaningful  only  if 
the  error  variance  is  less  than  tolerance 
level, 

Var(z)<^  (12) 

representing  an  ellipsoid  S^f)  with  y(f)  as  its 
center  (Fig.  1).  VPP  is  represented  by  a 
time  period  (f  -  f0)  at  which  z  (the  error)  first 
goes  out  of  the  ellipsoid  S^f).  The 
parameter  (f  -  f0)  is  usually  a  random 
variable.  The  joint  probability  density 
function  of  (f  -  t0)  and  z0  satisfies  the 
backward  Fokker-Planck  equation  [Section 
3.6  in  Gardiner,  1983] 


2 


dp 

d~t 


[f(z0,0_ 


dp 
d  zn 


1  2  2  d2P 

—q  z  - 

2  0  3  z  3  z 


=  0. 


(13) 

If  the  initial  error  exceeds  the  tolerance  level 
[i.e.,  z0  hits  the  boundary  of  Se  (to)],  the 
model  loses  prediction  capability  initially 

P((0,zoeS£((0),0)  =  0.  (14) 


The  model  also  loses  prediction  capability  at 
t  =  °°  [Gardiner,  1984] 

Lim  P(to,zQ,t-to)  =  0  . 

t— »°° 

Temporal  integration  of  P(t0,  z0,  t-  t0)  from  f0 
to  °°  should  be  one, 

]p(t0,z0,f-t0)dt  =  l.  (15) 

Since  P(t0,  z0,  t  -  t0)  is  the  probability  of 
VPP  (f  -  t0)  with  the  initial  error  vector  z0  at 
t0,  its  first  moment 


Ti(z0)  =  \P(tO’Zo’t-to){t~to)dt 

tn 

(16) 

‘0 

denotes  the  ensemble  mean 

VPP.  Its 

second  moment 

T2(zo^  =  ]P^tO’zo’t-to)(t-to)  dt 

(17)  indicates  the  variation  of  VPP.  Both 
(16)  and  (17)  indicate  the  dependence  of  the 
ensemble  mean  and  variance  of  VPP  on  the 
initial  error  z0. 


X 


Figure  1.  Phase  space  trajectories  of  model 
prediction  y  (solid  curve)  and  reality  x 
(dashed  curve)  and  error  ellipsoid  S^t) 
centered  at  y.  The  positions  of  reality  and 
prediction  trajectories  at  time  instance  are 


denoted  by  “*”  and  “o  ”,  respectively).  A 
valid  prediction  is  represented  by  a  time 
period  (t  -  to)  at  which  the  error  first  goes  out 
of  the  ellipsoid  S^t). 


3.3.  Autonomous  Dynamical 
System 

Usually,  two  steps  are  used  in 
computing  the  mean  and  variance  of  VPP: 

(a)  to  obtain  P(t0,  z0,  t  -  t0)  after  solving  the 
backward  Fokker-Planck  equation  (13),  and 

(b)  to  compute  the  mean  and  variance  of 
VPP  using  (16)  and  (17).  For  an 
autonomous  dynamical  system, 


we  multiply  the  backward  Fokker-Planck 
equation  (13)  is  by  (f  -  f0)  and  (f  -  f0)2  , 
integrated  both  equations  with  respect  to  t 
from  t0  to  °°,  and  obtain  the  mean  VPP 
equation 


...  .  3z\ 

f(zo)v-L  + 
ozn 


2  2  -j2, 

Q  Z0  O  T  _ 


2  dz0dz0 


=  -l 


(18) 


and  the  VPP  variability  equation 

3r,  ~2~2 


f(zo)v^  + 
oz„ 


C1  z  o  d  T2 

2  3z03z0 


=  -2r  (19) 


Both  (18)  and  (19)  are  linear,  time- 
independent,  and  second-order  differential 
equations  with  he  initial  error  z0  as  the  only 
independent  variable.  To  solve  these  two 
equations,  two  boundary  conditions  for  x-i 
and  x2  are  needed.  When  the  initial  error 
reaches  the  tolerance  level,  Var(z0)=  i,  the 
model  prediction  capability  is  lost  no  matter 
how  good  the  model  is,  and  X!  and  x2 
become  zero, 

r  =  0,  rn  =  0  for  Var^o)^.  (20) 


When  the  initial  error  is  below  the  noise 
level,  £n0ise,  the  initial  condition  is  considered 
as  ‘accurate’.  The  model  capability  of 
prediction  does  not  depend  on  the  initial 
condition  error  (i.e.,  X!  and  x2  are 
independent  on  z0), 


=  0,  ^  =  0  for  Var(z0)=£2n 

uZn  OZn 


(21) 


4.  Example 

4.1  One-Dimensional  Stochastic 
Dynamical  System 

We  use  a  one-dimensional 
probabilistic  error  growth  model  [Nicolis, 
1992] 


3 


^  =  (^-gf)  +  v(r)#,  0<^<oo  (22) 

dt 

as  an  example  to  illustrate  the  process  of 
computing  mean  VPP  and  VPP  variability. 
Here,  the  variable  £  corresponds  to  the 
positive  Lyapunov  exponent  a,  g  is  a  non¬ 
negative  parameter  whose  properties 
depend  on  the  underlying  attractor,  and  v(t )£ 
is  the  stochastic  forcing  satisfying  the 
condition 

(v(0}  =  0,  (v(t)v(t  '))  =  cfS(t-t'). 

Without  the  stochastic  forcing,  v(t)E  the 
model  (22)  becomes  the  projection  of  the 
Lorenz  attractor  onto  the  unstable  manifold. 

4.2.  IE  Analysis 

The  PDF  of  the  random  variable  £,  in 

(22) ,  P  f),  satisfies  the  Fokker-Planck 
equation  [Nicolis,  1992] 

df  +  -^[(^-gZ2)Pl  =  ^^(Z2P)  (23) 

at  oq  1  dq~ 

Multiplying  both  sides  of  (23)  by  £  and 

averaged  over  P  (£  f)  leads  to  the  time 
evolution  of  the  ensemble  mean  error 

^-  =  <r{$)-g{$)  -g{A£)  (24) 

2 

where  (A£  )  is  the  variance.  Both  (23)  and 

(24)  are  used  to  evaluate  model 
predictability  with  given  initial  condition, 

(£)  |t=o,  or  initial  error  distribution,  P  (£0). 

For  example,  Nicolis  [1992]  investigated  the 
predictability  of  the  stochastic  dynamical 
system  (5).  With  a  given  initial  error 
distribution 

P(£0)  =  *(£),  (25) 

she  integrated  the  Fokker-Planck  equation 

(23)  to  obtain  the  time  evolution  of  P(g,t) . 
With  a  given  initial  error 

(4)  |, =o=1 .6x1 0"2  (26) 

she  integrated  the  nonlinear  equation  (24) 
numerically  to  obtain  the  time  evolution  of 
ensemble  mean  error  {£,)  t.  The  model 
predictability  can  be  easily  evaluated  from 
temporal  variability  of  P  (£,  t)  and  (£)  t. 

4.3.  Equations  for  Mean  and 
Variance  of  VPP 

How  long  is  the  model  (22)  valid 
since  being  integrated  from  the  initial  state? 


Or  what  are  the  mean  and  variance  of  VPP 
of  (22)?  To  answer  these  questions,  we 
should  first  find  the  equations  depicting  VPP 
of  (22).  Applying  the  theory  described  in 
Sections  3.2  and  3.3  to  the  model  (22),  the 
PDF  for  the  random  variable  (f  -  f0),  [i.e., 
VPP]  satisfies  the  backward  Fokker-Planck 
equation, 


dP  r  ,  r2,dP  1  2d-P 

log  gg  L  £  .  q  e ,  =  0 

dt  0  0  df0  2  d^- 

(27) 

with  the  initial  error  (£0)  bounded  by, 

£  <E  <e 

7  noise  ~~  t>0  —  °  ■ 

Furthermore,  equations  (18)  and 
become  ordinary  differential  equations 

(19) 

( ^0  d  T 
(°5o  gg o)  0  ~TJ+=  1 

dg o  2  dq0 

(28) 

(^0-gO^p  +  ^f^  =  -  2r> 

dg o  2  dq0 

(29) 

with  the  boundary  conditions, 

Tx  =  0,  r7  =  0  for£0=£ 

(30) 

^  =  0,  4r  =  0for*=^noise. 

(31) 

4.4.  Analytical  Solutions 

Analytical  solutions  of  (28)  and  (29) 
with  the  boundary  conditions  (30)  and  (31) 
are 


T^^nois,£)  =  ^\y  ql  exp(^) 

q  So  q 

y  2<T_7 

f  xq  exp(- — f-x)dx  I dy 


(32) 


and 


1  _-o 


p<%o 

1  So  1 


y  ~x~2 

J  Tx(x)xq  exp(- — Y~x)dx 

f  .  CP 


dy 


(33) 


where 


£  =  £  Is ,  £  =£  .  le 

*  0  ^0  ’  ^  noise  ^  noise 


are  non-dimensional  initial  condition  error 
and  noise  level  scaled  by  the  tolerance  level 
e ,  respectively.  For  given  tolerance  and 
noise  levels  (or  user  input),  the  mean  and 


4 


variance  of  VPP  can  be  calculated  using 
(32)  and  (33). 

4.5.  Dependence  of  xi  and  x2  on 

o’  £  noise’  ^ ) 

To  investigate  the  sensitivity  of  r-i  and 
t2  to  B  n,B  ,  and  e,  we  use  the  same 

values  for  the  parameters  in  the  stochastic 
dynamical  system  (5)  as  in  Nicolis  [1992] 
o=0.64,  g=0.3,  <72=0.2.  (34) 

Figures  2  and  3  show  the  contour  plots  of 

.  ,£)  and  .  ,£)  versus 

0 7  ^  noise7  7  zw  07  •’noise7  7 

(B  n,B  )  for  four  different  values  of  e 

(0.01,  0.1,  1,  and  2).  Following  features  can 
be  obtained:  (a)  For  given  values  of 

(B  n,B  )  [i.e.,  the  same  location  in  the 

contour  plots],  both  x-,  and  x2  increase  with 
the  tolerance  level  e.  (b)  For  a  given  value  of 
tolerance  level  e,  both  x-i  and  x2  are  almost 

independent  on  the  noise  level  £„oise 
(contours  are  almost  paralleling  to  the 
horizontal  axis)  when  the  initial  error  (<^0 )  is 

much  larger  than  the  noise  level  (£noise)- 
This  indicates  that  the  effect  of  the  noise 
level  (B  )  on  xi  and  x2  becomes  evident 

only  when  the  initial  error  0)  is  close  to 
the  noise  level  ( ^ nniSt,  ).  (c)  For  given  values 
of  (e,  B,nuiS(, ),  both  x-i  and  x2  decrease  with 
increasing  initial  error  £ 0. 


Figures  2.  Contour  plots  of  ,£) 

versus  (E  ,E  )  for  four  different  values  of  e 

O’  r  noise' 

(0.01,  0.1,  1,  and  2)  using  Nicolis  model  with 
stochastic  forcing  q2  =  0.2.  The  contour  plot 
covers  the  half  domain  due  to  E  >E 

0  noise 

Figures  4  and  5  show  the  curve  plots  of 
T .  ,£)  and  .  ,£)  versus 

o7  j  noise7  7  07  •'noise7  7 

%  Q  for  four  different  values  of  tolerance 
level,  e  (0.01,  0.1,  1,  and  2)  and  four 
different  values  of  random  noise  B  .  (0.1, 

0.2,  0.4,  and  0.6).  Following  features  are 
obtained:  (a)  xi  and  x2  decrease  with 
increasing  ^  ,  which  implies  that  the  higher 

the  initial  error,  the  lower  the  predictability 
(or  VPP)  is;  (b)  x-,  and  x2  decrease  with 

increasing  noise  level  B  .  ,  which  implies 
that  the  higher  the  noise  level,  the  lower  the 
predictability  (or  VPP)  is;  and  (c)  xi  and  x2 
increase  with  the  increasing  e,  which  implies 
that  the  higher  the  tolerance  level,  the  longer 
the  VPP  is.  It  is  noticed  that  the  results 
presented  in  this  subsection  is  for  a  given 
value  of  stochastic  forcing  (q2  =  0.2)  only. 


Figures  3.  Contour  plots  of  t2(£  ,£) 

versus  (E  ,E  )  for  four  different  values  of  e 

o’  r  noise ' 

(0.01,  0.1,  1,  and  2)  using  Nicolis  model  with 
stochastic  forcing  q2  =  0.2.  The  contour  plot 
covers  the  half  domain  due  to  E  >E 

0  noise 
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4.6.  Dependence  of  xi  and  x2  on 
Stochastic  Forcing  (q2) 

To  investigate  the  sensitivity  of  r-i  and 
t2  to  the  strength  of  the  stochastic  forcing, 
q2,  we  use  the  same  values  for  the 
parameters  (a  -  0.64,  g  =  0.3)  in  (34)  as  in 
Nicolis  [1992]  except  q  ,  which  takes  values 
of  0.1, 0.2,  and  0.4. 

Figures  6  and  7  show  the  curve 

plots  Of  ()>?  noise’ q2)  3nd 

T2^^noise^2)  versus  <f0  for  two 
tolerance  levels  (e  =  0.1,  1),  two  noise  levels 
(E  .  =0.1,  0.6),  and  three  different  values 

of  q2  (0.1,  0.2,  and  0.4)  representing  weak, 
normal,  and  strong  stochastic  forcing.  Two 
regimes  are  found:  (a)  Xi  and  x2  decrease 
with  increasing  q2  for  large  noise  level 

( ^ noise =  0-6).  (b)  'Ci  and  x2  increase  with 

increasing  q2  for  small  noise  level  (E  .  - 

0.1),  and  (c)  both  relationships  (increase 
and  decrease  of  x-i  and  x2  with  increasing  q2) 
are  independent  on  e.  This  indicates  the 
existence  of  stabilizing  and  destabilizing 
regimes  of  the  dynamical  system  depending 
on  stochastic  forcing.  For  a  small  (large) 
noise  level,  the  stochastic  forcing  stabilizes 
(destabilizes)  the  dynamical  system  and 
extends  (shortens)  the  mean  VPP. 


0.6  0.8 


0.4  0.6  0.8  1 


Figures  4.  Dependence  of  zjj;  ,£)  on 

the  initial  condition  error  t  for  four  different 

5  0 

values  of  c  (0.01,  0.1,  1,  and  2)  and  four 


different  values  of  random  noise  E  (0.1, 

”  noise  '  ’ 

0.2,  0.4,  and  0.6)  using  Nicolis  model  with 
stochastic  forcing  q2  =  0.2. 


The  two  regimes  can  be  identified 
analytically  for  small  tolerance  level 

( £  — >  0 ).  The  initial  error  E,  should  also  be 
small  (E 0  ' £)■  The  solutions  (32)  becomes 


Lim  r  (£ 

£->0 


noise 


,£) 


1 

<7-q2/ 2 


2  a 

, ,2  - 1 

q  p  qi 

o  —  „2  noise 

2(7- q 


i 


^0 


-1 


(35) 


The  Lyapunov  exponent  is  identified  as  (er- 
q2/2)  for  dynamical  system  (22)  [Hasmin’skii 

1980].  For  a  small  noise  level  (Enoise'  1), 


the  second  term  in  the  bracket  of  the 
righthand  of  (35) 


R=- 


2(7 

n2  - 1 

q  p  ?2 

O  —  „2  noise 

2(7 -q 


f  i 


-1 


(36) 


is  negligible.  The  solution  (35)  becomes 


Lim  Ttfo,{n 


,£)  = 


1 


< 7-q 2/ 2 


-In 


1 

7^ 


(37) 


which  shows  that  the  stochastic  forcing  (q  ± 
0),  reduces  the  Lyapunov  exponent  (er  - 
q2l2),  stabilizes  the  dynamical  system  (22), 
and  in  turn  increases  the  mean  VPP.  On  the 

other  hand,  the  initial  error  E  reduces  the 
mean  VPP. 

For  a  large  noise  level  E  .  ,  the 

*  noise  1 

second  term  in  the  bracket  of  the  righthand 
of  (35)  is  not  negligible.  For  a  positive 
Lyapunov  exponent,  2 (7-q2  >0,  this  term  is 
always  negative  [see  (36)].  The  absolute 
value  of  R  increases  with  increasing  q2 

(remember  that  E noise  <1 ,  7o<  ”0-  Thus,  the 
term  (R)  destabilizes  the  one-dimensional 
stochastic  dynamical  system  (22),  and 
reduces  the  mean  VPP. 


6 


Figures  5.  Dependence  of  t2(E 0,£  n.  ,e)  on 

the  initial  condition  error  t  for  four  different 

=  0 

values  of  e  (0.01,  0.1,  1,  and  2)  and  four 

different  values  of  random  noise  E  (0.1, 

0.2,  0.4,  and  0.6)  using  Nicolis  model  with 
stochastic  forcing  q2  =  0.2. 


Figures  6.  Dependence  of  0,£ noise,q2) 

on  the  initial  condition  error  t  for  three 

r  0 

different  values  of  the  stochastic  forcing  qr2 
(0.1,  0.2,  and  0.4)  using  Nicolis  model  with 
Two  different  values  of  f  (0.1,  and  1)  and  two 

different  values  of  noise  level  E  (0.1,  and 
0.6). 


0.2  0.4  0.6  0.8  1  0.6  0.7  0.8  0.9  1 


Figures  7.  Dependence  of  T2(£ (),£ ^cf  ) 

on  the  initial  condition  error  E,  for  three 

different  values  of  the  stochastic  forcing  qr2 
(0.1,  0.2,  and  0.4)  using  Nicolis  model  with 
Two  different  values  of  e  (0.1,  and  1)  and  two 

different  values  of  noise  level  E  (0.1,  and 
0.6). 


6.  Conclusions 

(1)  The  model  valid  prediction 
period  (t  -  f0)  depends  not  only  on  the 
instantaneous  error  growth,  but  also  on  the 
noise  level,  the  tolerance  level,  and  the 
initial  error.  A  theoretical  framework  was 
developed  in  this  study  to  determine  the 
mean  (xi)  and  variability  (x2)  of  model  valid 
prediction  period  for  nonlinear  stochastic 
dynamical  system.  The  joint  probability 
density  function  of  the  valid  prediction  period 
and  initial  error  satisfies  the  backward 
Fokker-Planck  equation  when  the  valid 
prediction  period  is  assumed  homogeneous. 
After  solving  the  backward  Fokker-Planck 
equation,  it  is  easy  to  obtain  the  ensemble 
mean  and  variance  of  the  model  valid 
prediction  period. 

(2)  Uncertainty  in  ocean  (or 
atmospheric)  models  are  caused  by 
measurement  errors  (initial  and/or  boundary 
condition  errors),  model  discretization,  and 
uncertain  model  parameters.  This  leads  to 
the  inclusion  of  stochastic  forcing  in  ocean 
(atmospheric)  models.  The  backward 
Fokker-Planck  equation  can  be  used  for 
evaluation  of  ocean  (or  atmospheric)  model 
predictability  through  calculating  the  mean 
model  valid  prediction  period. 
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(3)  For  an  autonomous  dynamical 
system,  time-independent  second-order 
linear  differential  equations  are  derived  for 
Ti  and  t2  with  given  boundary  conditions. 
This  is  a  well-posed  problem  and  the 
solutions  are  easily  obtained. 

(4)  For  the  Nicolis  [1992]  model,  the 
second-order  ordinary  differential  equations 
of  Ti  and  t2  have  analytical  solutions,  which 
clearly  show  the  following  features:  (a) 
decrease  of  Ti  and  r2  with  increasing  initial 
condition  error  (or  with  increasing  random 
noise),  (b)  increase  of  p  and  r2  with 
increasing  tolerance  level  e. 
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